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First order phase transitions in the early universe can give rise to a stochastic background of 
gravitational waves. A hypothetical first order electroweak phase transition is particularly in- 
teresting in this respect, since the signal is in the good frequency range to be detectable by the 
space interferometer LISA. Three main processes lead to the production of the gravitational 

— h wave signal: the collision of the broken phase bubbles, the magnetohydrodynamical turbu- 

lence in the plasma stirred by the bubble collisions, and the magnetic fields amplified by the 
magnetohydrodynamical turbulence. The main features of the gravitational wave spectrum, 
such as the peak frequency, the amplitude, and the slopes both at low and high wave-number 

,l_> , can be predicted by general arguments based on the characteristics of the source: in particu- 

lar, the structure of its space and time correlation. We find that the gravitational wave signal 
from a first order phase transition occurring at electroweak symmetry breaking falls into the 
LISA sensitivity range if the phase transition lasts for about one hundredth of the Hubble 
time and the energy density of the turbulent motions is about twenty percent of the total 

^ ' energy density in the universe at the phase transition time. 
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1 Introduction 



—J . It is likely that in the years to come the interferometers LIGO^ and VIRGO ^ will provide 

the first direct detection of gravitational waves (GWs). Terrestrial interferometers operate in 

the frequency range 10-1000 Hz, while the space-based interferometer LISA will have its best 

sensitivity around the milliHertz frequency. A possible target of these experiments is a stochastic 

background of GWs of cosmological origin . The detection of such a background, relic of the 

early universe, would have a profound impact both on cosmology and on high energy physics. 

Once GWs have been emitted, in fact, they propagate freely through the universe, carrying 

direct information on the physical process that generated them: their detection would therefore 

provide a new way of probing the primordial universe, and correspondingly physics at very high 

energies, which would not be accessible otherwise. 

One of the processes that can generate such a stochastic background of GWs is a relativistic 

first order phase transition a l " ! ' l " l " l. In the course of its adiabatic expansion, the universe might 

have undergone several phase transitions driven by the temperature decrease. The nature of the 

phase transitions depends on the particle theory model, but if they are first order they proceed 

through the nucleation of broke n phase bubbles, which is a very violent and inhomogeneous 
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process capable of sourcing GWs U | ' '. In the following we review the different mechanisms 

of generation of GWs by a first order phase transition, and show that the characteristic fre- 
quency and amplitude of the GW signal are related respectively to the temperature at which 
the phase transition occurs and to its strength. We also review the analytical method developed 
in RefsJ ' ' \ which predicts the main features of the GW spectrum starting from a descrip- 
tion of the source which we tried to make as simple and as model independent as possible. The 



GW signal depends explicitly on a few free parameters: under s pecific ch oices for these parame- 
ters, it can be potentially interesting for observation with LISA ^ ' ' ' or, more speculatively, 
advanced LIGO. 

2 Gravitational wave background of cosmological origin 

In the cosmological context GWs are small perturbations of the FRW metric represented by the 
transverse-traceless tensor hij which is first order in cosmological perturbation theory: 

ds 2 = a 2 (t)(dt 2 - (5ij + 2h ij )dx i dx j ) , K = h^y = , (1) 

where we assume flat spatial sections and t denotes conformal time. Inserting the perturbed 
metric into Einstein's equations G^ u = 8ttGT^ u one obtains the evolution equation for the tensor 
perturbation, which in a radiation dominated universe takes the form 

KjiK t) + ^(k, t) + k 2 h i3 (k, t) = 8irGa 2 T^ T \k, t) , (2) 

where T*>- ' (k, t) is the transverse-traceless part of the energy momentum tensor sourcing the 
GWs, i.e. the tensor anisotropic stress. A source of GWs operating in the primordial universe 
is described as a stochastic process, and generates a stochastic background of GWs statistically 
homogeneous, isotropic, unpolarised and assumed to be gaussian. The energy density of the 
GWs, normalised to the critical energy density of the universe today p c , is given by the integral 
over wave-number of the power spectrum 

(hjjhv) _ fdk dQ GW d^ GW _ k 3 \h\ 2 

UGW ~ 87rGp c a2 ~ J k dlogfc' 6re d\ogk " 2(27r)3G Pc a 2 ' W 

and (h ij (k,t)h* j (q,t)) = (27r) 3 5(k-c l )\h(k,t)\ 2 . (4) 

If we consider a source of GWs active at a given time £* while the universe is in a phase of stan- 
dard FRW expansion (i.e. anytime besides inflation), then causality constrains the characteristic 
frequency of the emitted GWs to be larger than the causal horizon of the universe at the time of 
generation: k* > H*, where H* denotes the conformal Hubble parameter and fo* is the comov- 
ing wave-number. For example, GWs generated at the electroweak phase transition (EWPT) 
at T* ~ 100 GeV must have a characteristic frequency k* > 10~ 5 Hz, while the characteristic 
frequency of GWs generated at the QCD phase transition at T* ~ 100 MeV can be much lower, 
k* > 10~ 8 Hz. This estimation, based on the causality argument, shows that GWs generated 
during the EWPT are potentially interesting for detection with the space interferometer LISA, 
which operates in the frequency window from 10 -4 to 1 Hz. On the other hand, GW production 
at the QCDPT can fall into the frequency range of detection with pulsar timing array^^. We 
now proceed to analyse which kind of processes can act as sources of GWs during a relativistic 
phase transition occurring in the early universe. 

3 Gravitational waves from phase transitions 

There are a variety of processes related to primordial phase transitions that can lead to the 
production of GWs, as for example cosmic strings'^ or scalar field relaxation ^4 Here we 
concentrate specifically on the processes driven by bubble nucleation during a first order phase 
transition. The EWPT in the standard model is a crossover, and it is not expected to lead to 
any appreciable cosmological signal; however, deviations from the standard model in the Higgs 
sector or introducing supersymmetry can lead to a first order EWPT^^. Similarly, the QCDPT 



which is also predicted to be a crossover by lattice simulations'^, can become first order if 
the neutrino chemical potential is sufficiently large^l GWs detection would help to probe the 
nature of these phase transitions, and provide interesting information on the underlying particle 
theory. 

A first order phase transition proceeds through the nucleation of true- vacuum bubbles, which 
in the end of the transition collide and convert the entire universe to the broken phase. The 
collisions break the spherical symme try of the b ubble walls, generating a non-zero anisotropic 
stress which acts as a source of GWs' ' ' ' '. Moreover, bubble collision causes an injection 
of energy in the primordial plasma, which has a very high Reynolds number (of the order of 10 13 
at 100 GeV and at the typical scale of the bubbles'"^ : the energy injection leads to the formation 
of magnetohydrodynamic (MHD) turbulence in the fluid, which also sources GWs through the 
anisotropic stresses of the chaotic fluid motions . MHD turbulence also 

leads to the amplification of small magnetic fields generated by charge s eparation at the bubble 
wall, which have non-zero anisotropic stress and act as a source of G Ws l 25 | 26 | 27 | 28 | 29 | .S() ]3TITa. 

There are in summary three processes which can lead to the production of GWs towards the 
end of a first order phase transition. They are all related to the collision of bubbles, which in- 
volves two quantities: the duration of the phase transition, commonly denoted by the parameter 
P , and the typical size of the bubbles at the moment of collision, R* ~ v^ (3 , where Vb is the 
bubble wall velocity. The characteristic frequency of the GWs generated by the three processes 
can correspond either to the duration of the phase transition or to the bubble size: k* ~ (3 , 
R~ l , depending on the details of the time evolution of the source (c/. section d|). Assuming for 
the moment fc* ~ /3, one obtains the following order of magnitude estimate of the characteristic 
frequency in Hz: 

h =* 10~ 2 — * mHz . (5) 

H* 100 GeV V ; 

The parameter (3/%* represents the ratio of the duration of the phase transition with respect 
to the Hubble time. Since the entire universe must be converted to the broken phase, the phase 
transition must complete faster than Hubble time: a typical value for the EWPT isr /3/%* — 100. 
From Eq. ([5]) one gets that the characteristic frequency of GWs emitted at the EWPT is of the 
order of the milliHertz, and falls in the frequency range of the space interferometer LISA, which 
covers the frequency region 10 _4 Hz < k < 1Hz. LISA reaches its best sensitivity precisely 
around a milliHertz, where it can detect GWs with amplitude corresponding to fi^w — 10~ 12 . 
Starting from the perturbed Einstein equations SG^u = 8ttG5T^ u , one can give a simple 
order of magnitude estimate of the GW amplitude, which shows how the result depends on the 
duration and the energy density of the GW source. Since we are merely interested in determining 
the scaling of the GW amplitude, we rewrite Einstein equations simply as j3 2 h ~ 8ttGT, where 
h denotes the amplitude of the tensor perturbation, T the energy momentum tensor, and we 
inserted 1//3 as the characteristic time on which the perturbation is evolving (we have dropped 
indices for simplicity). This suggests that h ~ 8irGT//3, and the GW energy density becomes 
then /?gw ~ h 2 /(8irG) ~ 8ttGT 2 / ' j3 2 . For the three processes under analysis here, the energy 
momentum tensor of the GW source can be rewritten in all generality as T/p* ~ O*, where p* 
denotes the energy density in the universe at the GW generation time, and the parameter f2* 
denotes the relative energy density available in the source for the GW generation. We can now 
rearrange the equation for pew accounting for the fact that the phase transition takes place in 
the radiation-dominated universe, that the GW energy density also evolves like radiation, and 
inserting Friedmann equation 3H* = 8irGp* . It becomes then 
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^GW-^rad^yJ (SSr, (6) 

where f2 ra( j denotes the radiation energy density parameter today. The above equation shows 



that the GW energy density scales like the square of the ratio of the GW source duration and 
the Hubble time, and the square of the energy density in the source. Using /i 2 il ra d = 4.2 • 1CP 5 
and H*/f3 ~ 0.01, it turns out that a GW source with relative energy density of the order of 
O* ~ 0.1 could generate a GW signal above the lowest sensitivity of the space interferometer 
LISA. A detectable signal can therefore arise only from very energetic processes, which must 
involve at least 10% of the total energy density in the universe: namely, the phase transition 
must be strongly first order. 

4 The gravitational wave power spectrum 

Having demonstrated that GWs generated by a first order phase transition have a characteristic 
frequency and amplitude which can be interesting for detection, we now proceed to evaluate the 
main features of the GW power spectrum. The aim of this section is to show how the form of 
the power spectrum depends on the characteristics of the stochastic source, in parti cular the 
structure of its spatial and temporal correlation. It is based on the work done in Refs J ' ' ' . 
As derived in^, Eq. (|2|) can be solved in the radiation era (neglecting changes in the number 
of effective relativistic degrees of freedom), and once combined with definitions (J3J) and (|4|) it 
leads to the following formula for the GW power spectrum today : 

d^GW 4Q rad 3 /•*«" dh /•'«« dt 2 
where Tl(k,ti,t2) is the unequal time power spectrum of the source anisotropic stress, 

(n 4i (k,ti)n?.(q,t 2 )) = (2^) 3 <5(k-q)n(Mi,t2), (8) 

and we denote 11^ the dimensionless anisotropic stress: T^- ' = (p + p) Ily (p and p are the 

energy density and pressure of the primordial relativistic fluid, cf. Eq. (|2|) and Ref.LL^. The 
transverse traceless part of the spatial energy momentum tensor of the source contributes to 
the tensor anisotropic stress. In the bubble collisions case the spatial energy momentum tensor 
is given by Tjj = (p + p) ViVj/(l — v 2 ), where Vi is the velocity of the relativistic fluid at the 
bubble wall position^'. In the case of MHD turbulence, the energy momentum tensor can 
be decomposed into a part representing the turbulent velocity field and a part representing 
the magnetic field (for details, see^|). The turbulent velocity field part takes the form T^j = 
(p+p) ViVj (where we omit the gamma factor since we assume only mildly relativistic velocities), 
and the magnetic field one takes the form Tij = BiBj/^ir). We see that for all the three GW 
generation processes the tensor anisotropic stress is radiation-like, since the phase transition 
takes place in the radiation dominated universe. The parameter Q* introduced in section [3] is 
therefore given simply by f2* ~ (v 2 ) in the case of the turbulence and $7* ~ (b 2 ) in the case 
of the MHD processed magnetic field, where bi = B{/ \/~p^d is the dimensionless magnetic field 
parameter UM In the bubble case, the definition of fi* is a bit more involved, depending on 
whether the phase transition happens in a thermal bath or in vacuum: it is anyway related 
to the ratio of the kinetic energy density due to the bubble wall motions and the total energy 
density in the universe^'. 

Eq. ([7]) shows that the GW power spectrum is given by the double integral of the Green 
function of Eq. ([2]) multiplied by the unequal time anisotropic stress power spectrum, and 
its general shape as a function of wave-number k depends on the interplay among these two 
quantities. In particular, the spatial and temporal correlation of the source, together with 
its overall time evolution, determine the k and time dependence of the anisotropic stress power 
spectrum both at equal and unequal time, and in turns the GW spectrum. Much of the analytical 



work of RefsJ ' ' ' has dealt with the problem of modeling the statistical source for bubble 
collisions and MHD turbulence, and here we summarise the results obtained there. 

First of all, the /c-dependence of the equal time power spectrum Tl(k, ti,£i) is relatively easy 
to determine, following mainly from a causality argument. The bubble collisions and subsequent 
MHD turbulence are causal processes, characterised by a typical length-scale: the size of the 
bubbles at the moment of collision i?*, which also corresponds to the scale of energy injection in 
the primordial fluid generating the MHD turbulence. On scales larger than R*, the anisotropic 
stresses are not correlated and we expect the power spectrum to be flat, up to the wave-number 
corresponding to the inverse characteristic scale R* 1 . Beyond k ~ R*, the power spectrum 
decays with a slope that depends on the details of the source, and turns out to be A; -4 for bubble 
collisions'^, and A; -11 ' 3 for MHD turbulence!^. The spatial correlation structure of the source 
completely determines the /^-dependence of H(k, ti,ti), shown in Fig. Q] for the MHD turbulence 
case. 

Concerning the time dependence of H(k,ti,ti), it is due to the overall evolution of the GW 
source in time. In the case of bubble collisions, the source turns off right at the end of the phase 
transition, and lasts therefore for much less than one Hubble time. The overall time dependence 
of the source in this case comes mainly from the bubble expansion L^. In the case of MHD 
turbulence on the other hand, the dissipation of the turbulent motions in the primordial fluid 
is not very efficient due to the extremely low viscosity of the fluid itself!^. MHD turbulence 
is therefore maintained in the primordial fluid for several Hubble times, and the overall time 
dependence of the source depends both on the slow decay of the turbulent motions and on the 
growth of the turbulence characteristic scale associated with the decay. The influence of the 
total duration of the source on the shape of the GW spectrum shows up mainly at very large 
scales: if the source lasts several Hubble times, the GW signal is amplified on scales larger than 
the Hubble scale at the phase transition time k < 1-L *^. 

Moreover, in order to find the GW spectrum one needs to know the unequal time power 
spectrum of the anisotropic stress, as given in Eq. ([7]). This is much less obvious to determine, 
but one can identify a few simple models which are easy to dea l with an alytically and represent 
the main properties of the source. These have been discussed in! ' ' \ and in what follows we 
present the main results. In the bubble collision case, the source can be modeled as completely 
coherent, meaning that its time dependence is deterministic. In this case, the unequal time 
anisotropic stress power spectrum can be decomposed in terms of the equal time one as 



U(k,t u t 2 ) = ^U(k,t 1 ,t 1 )^U(k,t 2 ,t 2 ) . (9) 

This is a consequence of two main facts: first of all, the signal comes from the individual 
collision events, and each collision event can be assumed to be uncorrelated in time with the 
others; second, each collision event can also be assumed to be coherent, since the time evolution 
of the anisotropic stress related to the collisions is deterministic (it is basically only due to the 
growth in time of the bubbles'^). In the case of MHD turbulence, the situation is different. 
The turbulence can be viewed as a superposition of eddies of different size, each with its proper 
turnover time related to the eddy size. In the Kraichnan model!^, the eddy turnover time is the 
typical decorrelation time of turbulent motions. Therefore, the source has a finite decorrelation 
time which depends on the eddy size: this can be modeled with a top-hat decorrelation: 

n(fc, h,t 2 ) = {u(k, t 1: h)G[t 2 - h]&[i - k(t 2 - ti)] + h o t 2 } , (io) 

meaning that the source is correlated for time differences smaller than about one wavelength! 1 -^* 
\h-t 2 \ <k- 1 . 

The behaviour of the anisotropic stress unequal time power spectrum strongly influences 
the GW spectrum, and in particular the position of the peak. For a coherent source, the GW 
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Figure 1: Left panel: the anisotropic stress power spectrum at equal time as a function of the dimensionless 
variable kR* for the MHD turbulence. On scales larger than the characteristic scale R, the spectrum is flat 
since the source is uncorrelated, while on scales smaller than R* the spectrum decays as fc" 11 ' 3 (we assumed 
Kolmogorov turbulence). Right panel: the qualitative behaviour of the GW spectrum from the MHD turbulent 
source modeled as completely coherent (red, solid line) and as top-hat decorrelation (blue, dashed line). The 
peak position corresponds to the characteristic time-scale of the source fc* ~ /3 for the coherent case, and to 
the characteristic length-scale fc* ~ 1/i?* for the top-hat decorrelation case. The low frequency increase is 
k z by causality, and the high frequency slope depends both on the time decorrelation assumption and on the 

fc-dependence of the source power spectrum (left panel). 



spectrum becomes the square of the temporal Fourier transform of the source^', as can be seen 
by inserting Eq. ([9]) into ©. The source is characterised by the spatial correlation scale R* and 
the temporal correlation scale /3 _1 , related by R* = u&/3 . Since Vb < 1, one has R* < /3 _1 . On 
scales larger than both the characteristic spatial correlation scale i?* and temporal correlation 
scale /3 _1 , the space and time Fourier transforms of the source are constant because the source 
is not correlated (white noise). Therefore, for wave-numbers k < (3 < 1/R*, the GW spectrum 
simply increases as k 3 (cf. Eq. ([7])). However, for k > f3 the time Fourier transform is no longer 
constant and starts to decay as a power law, the exponent depending on the time differentiability 
properties of the source^^. In the bubble collision case, it turns out that this implies a k~ l 
decay for the GW spectrum at intermediate scales /3 < k < \/R*: the peak of the GW spectrum 
corresponds therefore to the characteristic time of the source^, k* ~ j3. For a source with finite 
decorrelation time as in Eq. (|10p . on the other hand, the GW spectrum is no longer related to 
the temporal Fourier transform of the source, and the situation changes. The spectrum still 
increases as k 3 on large scales, but the characteristic time of the source /3 _1 does not lead to 
any feature in the spectrum, which peaks instead at the inverse characteristic scale of the source 
fc* — R* 1 . The qualitative shape of the GW spectra coming from a coherent source and a 
top-hat decorrelating one is shown in Fig. [TJ 



5 Results 



Fig. [2] s hows the GW spectra generated by bubble collisions and MHD turbulence, as derived 
in RefsJ ' ' '. The parameter representing the duration of the phase transition is set to 
Tl*/ ' j3 = 0.01, and the parameter representing the relative energy density available in the source 
for the GW generation is set to O* = 0.2 in both cases (cf. Eq. [6]). This high value of O* implies 
a strongly first order phase transition, for which the vacuum energy density is about one third 
of the radiation one a = 1/3, and the bubble wall speed is close to the speed of light, Vb = 0.87. 
Correspondingly, the mean velocity of the turbulent motions must be of the order of the speed 
of sound, (v 2 ) = 1/3, and equipartition is assumed among the kinetic and the magnetic energies 
in the turbulence such that (b 2 ) ~ {v 2 ). 

In the two GW spectra of Fig. [2] we can distinguish the features presented in the previous 



sections. Both spectra rise as k 3 for small wave- numbers: this is simply the phase space volume 
(c/. Eq. ([7])) combined with the fc-dependence of the anisotropic stress power spectrum, which 
is flat at small wave-numbers because of causality. In the GW spectrum generated by bubble 
collisions the k s slope is maintained up to the peak, while in the one generated by MHD tur- 
bulence the slope changes to k 2 at sub- horizon scales k > 7i*: this difference is due to the fact 
that bubble collisions are a short lasting source, while MHD turbulence acts a source of GWs 
for several Hubble times^^. The GW spectrum due to bubble collisions peaks at a wave-number 
corresponding to the duration of the source fc* ~ j3, since the source is modeled as completely 
coherent in its time decorrelation; while the GW spectrum due to MHD turbulence peaks at a 
wave-number corresponding to the characteristic length-scale of the source, given by the bubble 
size at the end of the phase transition: k* ~ 2-k/R*, since the source is modeled following the 
top-hat Ansatz and has a finite decorrelation time corresponding to the eddy turnover time^^. 
At frequencies smaller than the peak, the GW spectrum from bubble collisions decays as A; -1 : 
this decay is related both to the fact that the source is coherent and to the thin wall approx- 
imation, which has been inserted in the analytical analysis to recover the result of numerical 
simulations ^ ' '. The GW spectrum from MHD turbulence decays at high frequencies with 
slopes that also depend both on the top-hat decorrelation structure and on the source power 
spectrum: they turn out to be k~ 5 ' 3 for the kinetic turbulence (coming directly from the as- 
sumption of a Kolmogorov spectrum) and k~ 3 ' 2 for the magnetic field (due to the assumption 
of an Iroshnikov-Kraichnan spectrum) L^. 

Fig. [3] shows the final GW power spectrum given by the sum of the bubble collisions signal 
and the MHD turbulence signal for the EWPT at T* ~ 100 GeV. The parameters are again 
~H*//3 = 0.01 and $7* = 0.2, so that the phase transition is assumed to be strongly first order. 
The signal is compared with the sensitivity curves of LISA and BBO: for a strongly first order 
EWPT the signal falls into the sensitivity range of LISA. Clearly this can only be achieved in 
the context of theories that go beyond the standard model of particle physics. Fig. [3] also shows 
an even more speculative case of a phase transition occurring at temperature T* = 5 • 10 6 GeV, 
with H*/(3 = 0.02 and SI* = 0.2: in this case, the signal could be interesting for advanced LIGO 
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Figure 3: Left panel: the total GW spectrum for bubble collisions and MHD turbulence generated at a strongly 
first order EWPT as a function of frequency in Hz. The duration of the phase transition and the energy density 
of the sources are given the values H*/fi — 0.01 and Q.^ = 0.2, and the signal falls into the sensitivity range of 
the space interferometer LISA. Right panel: the total GW spectrum from MHD turbulence only, generated at a 
strongly first order EWPT (left black curve) and at a hypothetical first order phase transition at temperature 
T* = 5 ■ 10 6 GeV. The parameters are H*/0 = 0.01 and Q* s = 0.2 for the EWPT and H*//3 = 0.02 and Q,* s = 0.2 

for the PT at T, = 5 ■ 10 6 GeV. Figure from 7 . 
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